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The rate equations of Statz and Dc Mars giving the time development of 
the inversion and photon number in a maser or laser are discussed ana- 
lytically with the aid of a mechanical analogy in which a particle moves in 
a potential well under tlie influence of a viscous damping force. The coor- 
dinate of this particle is analogous to the logarithm of the light output of the 
laser, and the amplitude, period, and damping of the motion can be directly 
related to the parameters of the rate equations. Simple analytic approxima- 
tions are developed for all of the quantities of experimental interest in the 
spiking pattern of a laser. Four relationships are given, which do not con- 
tain any of the rate equation parameters, whereby a spike pattern can be 
tested to determine if it is consistent with the usual rate equations. Systematic 
procedures are described for extracting all of the information contained in 
spike patterns. 

I. INTRODUCTION 

The most fruitful approach for the discussion of maser and laser 1 
behavior has been through rate equations describing the time rates of 
change of the atomic populations and the photon numbers of the electro- 
magnetic field. Bloembergen 2 introduced rate equations for the popula- 
tions in a paramagnetic maser and based his discussion on the steady- 
state solutions without explicitly considering the photon field. On the 
other hand, Shimoda, Takahasi, and Townes 3 have considered the pho- 
ton rate equations and on this basis have given a theory of maser ampli- 
fication without explicitly considering the atomic populations. Statz 
and De Mars 4 have shown that the transient behavior of masers depends 
upon coupled rate equations for both the populations and the photons. 
A number of authors have rederived these equations and discussed their 
applications to various maser systems. Considerable attention has been 
given to the question of whether these equations have periodic (un- 
damped) solutions. It has been shown by Makhov 5 and by Sinnett 6 that 
the small-signal solutions are always damped, and it has been pre- 
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suined and often confirmed by numerical computations that the same is 
true in the large-signal domain. Statz, et al. 7 have suggested that the 
damping of experimentally observed "spikes" in the output of lasers is 
probably sensitive to coherence and noise conditions, and not neces- 
sarily related to the damping predicted by the usual rate equations. 
They also suggest that the complicated spiking patterns frequently 
observed 8 are due to oscillation in many modes of the laser cavity. De- 
spite these doubts which have been cast on the adequacy of the Statz- 
De Mars rate equations for describing laser behavior, they still remain 
the logical starting point for any discussion of the power output of lasers, 
whether it be transient or steady-state, and of the dependence of power 
on the quality of the cavity, the intensity of the pumping light, or the 
linewidth, relaxation time, and concentration of the active atoms. 

Several ingenious suggestions have been made for modifying the rate 
equations so as to obtain periodic solutions. Statz and De Mars 4 and 
Makhov 5 propose that periodic solutions are obtained if terms are added 
to the rate equations representing cross relaxation in the inhomogene- 
ously broadened maser transition. On the other hand, Shimoda 9 sug- 
gests that periodic spiking can result if the losses due to absorption in the 
cavity can be partially saturated by the buildup of laser oscillations. Al- 
though these and other modifications may ultimately prove to be justi- 
fied and necessary in laser theory, we shall confine our attention in this 
paper to the original Statz-De Mars equations which relate the photons 
in a single cavity mode to a single quantity, the inversion, describing the 
atomic populations. We shall make it our task to understand as fully as 
possible the damped oscillatory solutions of these equations and how they 
may be applied to the study of the spiking phenomenon seen 8 in solid- 
state lasers. 

Rate equations in the simple form have been successfully applied by 
McClung and Hellwarth 10 and by Vuylstekc" to the giant-pulse laser, 
which produces a single very short and very intense burst of radiation. 
Wagner and Lengyel 12 have shown that an exact analytic solution can be 
obtained to a simplified rate equation which neglects spontaneous emis- 
sion and pumping during the pulse. By also neglecting the loss of pho- 
tons in the cavity Dunsmuir 13 has obtained a still simpler analytic solu- 
tion which is applicable to the rising portion of a spike or a giant pulse. 
We shall not consider further these exact solutions or the giant-pulse 
laser in this paper, but confine ourselves to those solutions representing 
repetitive pulsations, or spikes, in the ordinary laser. It is in this field 
that the greatest need now exists for an analytical discussion of the solu- 
tions of the rate equations. 
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A number of authors have endeavored to put the rate equations on a 
firmer theoretical basis. Following Anderson 14 and Clogston, 15 who first 
described masers in terms of the density matrix, treatments using the 
density matrix to derive rate equations have been given by Fain, et al., 16 
Kaplan and Zier, 17 and Pao. 18 By considering directly, without explicit 
use of the density matrix, the correlation functions of the electromag- 
netic field which are measured in simple maser experiments, McCumber 19 
has shown that the maser medium acts like a dielectric of negative con- 
ductivity; he concludes that the field and the dielectric satisfy rate equa- 
tions of the usual form, providing that the populations change by a small 
fractional amount during a coherence time (reciprocal linewidth) of the 
atomic system. Another formal theory, based on a successive approxima- 
tion approach to the quantum mechanical equations of motion, has been 
applied by Haken and Sauermann 20 to the frequency shifts and interac- 
tions of cavity modes in the laser. An extensive survey with bibliography 
of the early formal work has been given by Lamb. 21 In the present paper 
we shall not go into the theoretical basis for the rate equations, but con- 
fine ourselves entirely to the problem of solving the equations in the 
large-signal domain. 

Despite the fact that the rate equations are generally accepted, and 
that the general nature of the solutions has been familiar from digital 
computer calculations for some time, 5131722 it is still quite inconvenient 
in any particular case to compare observed spiking patterns in solid- 
state lasers with predictions of the rate equations. It has been neces- 
sary either to use the small-signal solutions and hope that they are not 
too inaccurate in the large-signal domain, or to resort to machine calcu- 
lations and try to fit three or more parameters to the data by trial and 
error. To be sure, much of the data on spiking is not amenable to analysis, 
consisting apparently of random spikes with widely varying amplitude, 
duration, and interval. Nevertheless, several laser systems are now known 
to give very regular pulsations of the type that might be consistent with 
the rate equations. Regular spiking patterns have been observed in 
OaF 2 :U +3 by Sorokin and Stevenson 23 and by O'Connor and Bostick, 24 
and in CaW0 4 :N<i +3 by Johnson and Nassau. 25 Recently the effect has 
also been seen in a highly perfect ruby by Nelson and Remeika, 26 and in 
a confocal ruby by Johnson, et al. 27 More extensive studies of highly 
regular spiking have been reported by Giirs 28 and by Hercher. 29 Thus it 
is clear that good data on spiking patterns can be obtained, and it 
therefore becomes cogent to inquire into practical and convenient means 
for analyzing this data and obtaining information from it. 

The information contained in spike patterns is of two distinct kinds, 
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which we may call qualitative and quantitative. Qualitatively, we can 
determine quickly by means of relationships given here whether a spike 
pattern is consistent with the rate equations. If we determine that cer- 
tain patterns are not consistent, we are spared the waste of time that 
would result from attempting to fit these patterns numerically by trial 
and error. It is not obvious at this writing that any of the regular pat- 
terns that have been reported are consistent, because the application of 
the consistency relationships requires a measurement of the ratio of the 
peaks to the valleys in the light output, and only the peaks are seen in 
pictures published so far. Thus we suggest that by extending spike studies 
to include the valleys new and interesting information can be obtained. 
It is to be expected that patterns which are regular but yet not consist- 
ent will turn out to be physically the most interesting of all, since they 
will point the way to new understanding of laser behavior. The quanti- 
tative information can only be obtained from patterns which are con- 
sistent, at least in some average sense, and consists in obtaining values 
for the physical constants which appear in the rate equations : 
N = the number of active laser atoms in the optical cavity 
t r = the relaxation time of the upper laser level, usually due to spon- 
taneous emission 

t p = the photon lifetime in the laser mode of the cavity 

t,„ = the mode time, the time for spontaneous emission into the laser 

mode 

t„ = the ground state time, the time spent by an atom in the ground 
state before being excited by the pump 

s = the source strength, the rate of production of laser photons by 
spontaneous emission, the pumping light, or any other noise source in 
the cavity, or any signal applied to the cavity. Although s may vary 
with time, it is convenient here to consider it with the constants. 

In principle all of these quantities except s could be directly measured 
or calculated from independent measurements on the laser material, the 
cavity, and the pump. The spiking data would then serve as confirmatory 
evidence. In many cases, however, spike patterns may prove to be the 
most convenient method of measurement. The last quantity, .s, is in 
some respects the most interesting. The first assumption would be that 
s is due entirely to spontaneous emission into the laser mode; if so, s 
could be calculated and comparison with the measured value would re- 
veal the verity of the assumption. Experiments could be carried out with 
an external weak signal from a monochromator to test the response of 
the laser as observed in its spiking patterns. Thus we hope that the 
analysis given here will help to stimulate new experiments by making 
rate equation analysis more convenient for the experimentalist. 
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II. FORMULATION OF THE RATE EQUATIONS 

The laser is a system consisting essentially of an optical cavity with 
very high Q in a few modes, low Q in all other modes, the laser medium 
containing the active atoms, and a pump, usually an intense light source, 
to excite the atoms into a broad hand of excited states. It is a property 
of the laser medium that the atoms decay from these states in an ex- 
tremely short time by nonradiative processes to one or more very sharp 
excited states, called the upper laser levels. The upper laser levels can 
decay radiativcly to a sharp lower level, called the lower laser level, 
which may be the ground state. If the lower laser level is not the ground 
state, we shall assume that very rapid nonradiative decay processes 
return the atom to the ground state. Thus we may always neglect the 
population of the pumping band and of the lower laser level if the latter 
is not the ground state. The laser transition takes place from the lowest 
of the upper laser levels, but it may sometimes be necessary to take into 
account the populations of nearby levels in thermal equilibrium with this 
level. The statistical weights of the laser levels must also be taken into 
account. 30 

The rate equations may always be written in the form (p = photon 
number, n = inversion) 

<lt tp tm 

dn Mo — m pn , x . 

-j- = — a,— (2) 

dt to t m ' 

as long as we consider only a single mode of the cavity, the laser mode, 
and neglect all atomic populations except the upper laser level and the 
ground state. The time t t , might be called the pumping relaxation time 

'•() tg t r 

since it represents the characteristic time in the response of the popula- 
tion inversion n to the pump. The inversion n may always be written 

n = N„ - (\/'w)N,, (4) 

where N u , Ni are the populations of the upper and lower laser levels 
respectively and w is the statistical weight of the lower relative to that 
of the upper laser level. In view of our assumptions, w will come in only 
when the lower laser level is the ground state. Let us suppose that the 
upper laser level is in thermal equilibrium " with certain other states 
not directly involved in the laser transition such that there is a tempera- 



1510 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1964 

ture-dependent probability P that an excited atom is in the upper laser 
level. Then we have for a three-level system, in which the lower laser level 
is the ground state, 



o = P + - 

w 



wo = PN t 



1 \ (5> 



and for & four-level system, in which the lower laser level is not the ground 
state, we have 

a = P 

(6) 
wo = PN t /t g . 

If the relaxation of the upper laser level is predominantly by spontane- 
ous emission to the lower laser level, there is a simple relation between 
t T , t m , and the linewidth Av (cps) of the laser transition 

t m = {%Trv l nJ/c 3 )Vt r Av, (7) 

where V is the volume of the cavity and w re f is the refractive index. This 
relation should not be taken too literally, since it takes no account of the 
anisotropy of the laser medium or the polarization properties of the 
transition, and V would have to be replaced by a suitable effective "op- 
tical volume" if the laser material does not fill the cavity. Nevertheless, 
it points out the important fact that t m varies with temperature in the 
same way as Ap and therefore is subject to control in spiking studies. 
Another constant subject to convenient control is t, , since \/t g is pro- 
portional to the pump intensity. Even when flash lamps are used for 
pumping it is still approximately valid to assume t is constant, since the 
time constant for the flash will usually be much longer than the interval 
between spikes. To assure that this is true, it would be advantageous to 
have the spike pattern commence when the flash is at its maximum in- 
tensity. It is not valid to assume without investigation that l/t„ is 
proportional to the total energy dissipated in the flash. Spiking data 
should always include a record of the flash as a function of time and an 
indication of when the spike pattern occurred. Also subject to experi- 
mental control is t p , the photon lifetime in the laser mode of the cavity. 
Presumably the losses in a good laser cavity can be estimated rather re- 
liably, so that t p can usually be directly calculated. The total number N 
of active laser atoms can ordinarily be determined in a given sample, but 
N does not appear to be a convenient parameter to vary in laser experi- 
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ments. The signal or source strength .s is best determined from an analysis 
of the spiking data as described in this paper. Although s = s(t) will in 
general be a function of time, it will be shown that the spiking pattern 
depends only on the value of 6- at the instant when the net losses in the 
laser mode vanish due to the buildup of inversion. We shall denote this 
time by U and the critical inversion by n x , where 

»1 = Im/tp (8) 

is the celebrated Schawlow-Townes 1 criterion for the buildup of laser 
oscillation. 

We now introduce dimensionless variables and parameters; in terms 
of the variables 

T = t/t p 

v = n{t P /t m ) = n/nj (9) 

P = p(a to/t m ), 
and the parameters 

w = tp/h 

£ = no(tp/tm) = n /rii (10) 

a = s(a tptu/t m ), 
the rate equations (1), (2) become 

p = _- = a — p -j- pq (11) 

(It 

v = p T = *(Z-v- pv). (12) 

Here p represents the photons, -q the inversion, and r the time, while w 
represents the pumping rate, £ the limiting inversion toward which the 
pump is tending to drive the system, and a the source. We see that there 
are really only three parameters in the rate equations; it follows that 
three relationships among the six relevant physical parameters with 
which we started can be obtained from spiking studies. Although we 
shall assume in our analysis that w and £ are constant, our results will 
provide a useful adiabatic approximation for the case of slowly varying 
a), £. Typical values of the parameters for a ruby laser will be given in 
the discussion of a numerical example. For the present we need only 
mention that a is relevant only in the initial growth of p prior to the on- 
set of laser gain. 
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If we regard a as a constant, the steady-state solution of ( 1 1 ) and (12) 

is 

Poo = }{(« + a - 1) + [(£ + « ~ I) 2 + 4<r]*} 
>/« = £{(£+*+ 1) " [(* + »- l) 2 + 4cr] J }. 

The sign of the radical is determined by the requirement that p ^ 0. 
In the limit a — * we obtain two cases, depending on whether £ < 1 or 
£> l.For£ < 1 

p M -> «r/(l - *) 
ff -"° £ < l; (14) 

r/oc -> f(l - p«) 

this is the case in which the limiting inversion n Q is less than n x given by 
(8), and laser oscillation does not occur. For £ > 1 

P. -» I - 1 

'"*° e>ij (is) 

17oo — > 1 

this describes the steady state of laser oscillation, which is usually ap- 
proached through a series of sharp pulses in p(r) called relaxation oscilla- 
tions, or spikes. 

III. THE NATURE OF THE SPIKING SOLUTIONS TO THE RATE EQUATIONS 

Iii this section we shall consider the nature of the solutions to (11) 
and (12) when £ > 1 and the initial conditions on p and ij correspond to 
very few photons and a small inversion n « Wi . This may be contrasted 
with the situation in the giant-pulse laser in which immediately after 
switching n » ni . We shall also assume the spiking condition 

«G«1, (16) 

which will be satisfied whenever spiking can be observed. It will be ap- 
parent later that when (16) is not satisfied there will be no spikes, but p 
will smoothly approach p M . This is the case in gas lasers, where the pump- 
ing rate \/t a has to be very high to overcome the high relaxation rate 
l/t r ~ 10 8 sec -1 . 

In view of (16) and (12), ij(t) will increase slowly with time and p 
in (11) can be neglected ; thus we have 

p ( T ) ~p( r ) ■ er(r)/(l - i(r)) (17) 

until r approaches n and t\ approaches unity. The behavior of p(r) is 
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shown in Fig. 1 for the initial condition p(0) = 0. Initially p rises with 
slope o-(O) and asymptotically approaches the adiabatic solution p(r) 
which it follows for a relatively long time until r — > ti . It follows that 
the initial condition on p(r) is unimportant. As r passes through ti , 
where 

r?(ri) = 1, (18) 

p(t) remains finite and is no longer given by (17); we then leave the 
adiabatic phase and enter the spiking phase of the time development of 
p(t). According to (11) the slope of p(r) at ti is <t(ti) = a x . Thus we 
construct an approximate solution by smoothly joining (17) to a line 
of slope o-i as shown in Fig. 1 . It follows that 

p ( Tl ) = pi = 2cn/(«0)», (19) 

where 



= e- i. 



(20) 



This is our first important result. It shows that the source occurs in 
spiking theory only as a parameter, the single value oi . Thus we may 
drop the subscript on a and let o\ = a, a, constant. 

Once we enter the spiking phase, 77(7-) remains close to unity, fluctuat- 




Fig. 1 — The adiabatic phase of the growth of p(r). The curve is the adiabatic 
solution (17), and the line segments of slope a are constructions to approximate 
p(t) near t = and t = ti where 9 = jji = 1. 
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ing up and down and eventually settling down to its steady-state value 
rj x = 1 . Thus we make the approximations 

rt - 1 « In 7,, (21) 



and 



t- l«/3+*(l - r,). (22) 

The rate equations (11), (12) can now be written 

— In p = In 7j (23) 

l-kiTj = »G8 - 8 In* - p). (24) 

or 

Since the inversion is not directly observed we eliminate In y from (23), 
(24) ; the result is most conveniently written 

* = w/3(l - c*) - wf* (25) 

in terms of the logarithmic light output 

* = In (p/j8). (20) 

The discussion of (25) is greatly facilitated by a mechanical analogy 
which is shown in Fig. 2. We regard ^ as the coordinate of a particle of 
unit mass moving in a one dimensional potential field. 

rw = -co/3 f d - ••) d* 

J o (27) 

= co/3(e* - * - 1). 

There is also a dissipative resistive force u$ as if the particle were mov- 
ing through a viscous medium. For the moment let us disregard the 
viscous force, in which case the total energy E of the particle is con- 
served 

E = V(V) + \¥. (28) 

The particle executes a periodic motion between extreme points >K, < 
and Wu > such that 

V(* m ) = V(*h) = E. (29) 

In the spiking phase it is permissible to neglect exp (^ m ). From (27) 
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Fig. 2 — The spiking phase of the time development of p(r) considered in terms 
of a mechanical analogy in which a particle with coordinate* = In (p//3) moves in 
a potential V(p) given by (27). The extremes of the motion *,„ and Vm are deter- 
mined by the total energy E of the particle. 

and (29) we then obtain a simple relation between ty m and ^^ 

&M - $m = exp (*jf). (30) 



This is our second important result. 
Let us define 



(31) 



*, = In (pi//3); 
then it follows from (28) that 

*« - *i - Ui/p,)72co0, (32) 

since ^i =» o"i/pi • It will be clear from the numerical example in the next 
section that the second term can be neglected, and we can write 

*n » *l • (33) 

This says that the small kinetic energy of the particle at the start of the 
spiking phase can be neglected. 

It is convenient to denote the successive times when tj(t) = 1 by 
i"i , t% , n , • • • . From (23) we see that these times correspond to ex- 
trema in the motion of ^; according to this convention the successive 
minima and maxima of ^ are 



minima: ( v l',), v l':i , *5 , *7 , 
maxima : % , ¥4 , % , • • • . 



(34) 
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We have placed Vi in parentheses, since it is a minimum only in the 
mechanical analogy and not in the observed light output. In the absence 
of damping, of course, we would have ^i = ^3 = ^6 = ■ ■ ■ , and 
>(r 2 = ^ 4 = ^ 6 = ■ • • . If (16) is satisfied, the damping will be small 
enough so that the motion is approximately periodic, or quasi-periodic. 
Thus to a good approximation Ave can compute the first maximum ^ 2 
from ^1 by means of (30). Before considering the damping we shall con- 
sider other quantities of experimental interest which are characteristic 
of the periodic motion. 

The maximum velocity $ mu = ^0 is capable of being measured experi- 
mentally from spiking patterns in which the spikes are well resolved in 
time. It follows from (28), (29) and (33) that 

*W = *o = [2co/3(e*' -*!- 1)] ! 

. (35) 

fc<P«tf<-*i- 1)] } . 

In numerical applications this can be used to ascertain the validity of 
(21), since according to (23) 

^0 = hi 7? max « T?,„ax ~ 1- (36) 

Iii general we can write near ^i 

if = {2[E - FC*)]) 1 
w [2u>/3(* - *i)] J ; 

thus the time dependence near a minimum is given by 

¥(r) « *i + |w0(t - t,) 2 . (38) 

Let us denote by m the full width in time of the minimum measured 
between points e times the minimum in light output; this is the same as 
the full width of ¥(t) measured between points V m + 1. Thus the dura- 
tion of the minima according to (38) is 

m = (8/co/3) J . (39) 

This applies to all minima regardless of damping, which provides a very 
convenient means for determining almost by inspection whether or not 
a spike pattern is consistent with the assumption of constant to/3. Even 
if the minima are not observed to have the same durations, it may be 
meaningful, in the sense that our theory provides an adiabatic approxi- 
mation, to apply (39) to each minimum separately and deduce the varia- 
tion of co/3. 



(37) 
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(40) 



Near the maximum ty 2 we write instead of (37) 

* « [2w/3(e* 2 - 1)(* 2 - ¥)]* 

» [2«/3(¥ 2 - *i - 1)(*2 -*)]*, 

where use has been made of (30); thus 

¥(t) « * 2 - |«|S(«| - *i - 1)(t - r 2 ) 2 . (41) 

The duration M of a maximum will be defined as the time interval meas- 
ured between points at l/e times the maximum in light output; this is 
the same as the interval between points at ~&m — 1. We see from the 
factor (^2 — ^i — 1) in (41) that .1/ will depend upon damping; it 
may be written in a general way 

M = >?)(*.„ - *„, - l)" 1 . (42) 

There is in this relation a certain ambiguity which is inherent in our 
method of regarding the motion as quasi-periodic. In applying the rela- 
tion we may wonder whether the minimum is the one preceding or fol- 
lowing the maximum. Within the accuracy of the quasi-periodic approxi- 
mation it makes no difference: either may be used, or the average of the 
two. 

The most readily observed quantity in spiking experiments is the in- 
terval between spikes, which can be identified with the period of the 
quasi-periodic motion 



I = <b (Nr/V, (43) 

where the integral is over one cycle. To evaluate / we use the approxima- 
tions (37) and (40), which are accurate near the turning points ^i and 
ty> respectively where 1,AJ> is large. Upon comparing (35) and (37) we see 
that (37) is reasonably accurate even at SP = providing | ^ | » 1. 
However, (40) is only accurate near ^ 2 , say in the range 

*2 - 1 < * ^ ^2 • 

Thus we shall use (37) in the range ty x ^ ^ ^ V,. and (40) in the range 
^< < ^ ^ ^2 , where ^,. is a crossover point which will be determined 
presently. The integral (43 ) can now be carried out to obtain 



/(¥ c ) « m 



(*■> - tf f )l 1 
'' •''■■ 7^-^-TijiJ. (44) 



Since both our approximations tend to underestimate 1/ty, we must 



(46) 
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choose ^ c so as to maximize /, which gives the condition 

Vc = *2 - 1. (45) 

Thus the interval I between spikes is given by 

/ = m [(p M - ■*„, - 1)* + (*m - ¥- - l) _i ] 

= M(* M - * m ). 

It is logical in this case to choose the minimum between the two maxima 
between which I is measured. There is still an ambiguity, however, in 
the choice of maxima. Since our approximation tends to underestimate 
7, it is good to use the larger of the two maxima. 

We now return to the equation of motion (25) and consider the damp- 
ing force —ca&r. If (16) holds, the damping will be small, and can be 
computed from the work done per cycle against the damping force. 

W = io£ (p 4fd9. (47) 



With damping present, energy is no longer conserved, but decreases by 
W every cycle of the motion until the particle eventually settles down at 
its equilibrium position ^ = 0, corresponding to the steady-state light 
output given by (15). Near ^i we can neglect e , so that (27) gives 

to/3(* 3 - *i) = W. (48) 

We evaluate (47) just as we did (43), using (37) and (40) with a cross- 
over point yp c , determined this time by the condition that W should be 
a minimum; the result is 

W = (4V2/3)(«e)(«0)*[(*J« - *,„ - 1)1+ (^ - * M - 1)*]. (49) 

Let us denote damping by the increment A^ m between successive minima, 
or A\F W between successive maxima. From (48) and (49) the general 
formula for damping of minima is 

M m = (4V2/3M(«0r*[(¥j, - **, - 1)» + (** - *™ " I)*]- (50) 

The choice of ¥ m is ambiguous, but V M refers to the maximum between 
the two minima of A^„, . It is obvious from the shape of the potential 
V(V) shown in Fig. 2 that the maxima will he less damped than the 
minima. For small damping we have 

u0(e* s - l)(¥ a -* 4 ) = W, (51) 

where W is now computed by integrating (47) from tf 2 around the cycle 
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and back to ^ 2 . The result is 

A*,, = (4V^/3)«€(«»8)" l [(*j, - *„, - 1)* + (*m - *«, - 1) _§ ] (52) 

with the familiar ambiguity in choice of ^ M . From (52), (46) and (39) 
we obtain the very convenient formula 

M„ = !(«£)/ (53) 

relating the damping directly to the interval. We note that all ambiguity 
lias disappeared from (53). 

We now have a complete arsenal of formulas with which to attack ex- 
perimental spiking patterns. Our formulas give all of the minima and 
maxima of ^ as well as the durations and intervals and the maximum of 
ty in terms of the dimensionlcss rate equation parameters cu, £, a and 
(3 = £ — 1. These formulas are valid in the spiking phase where 

¥* - *,„ - 1 » 1. (54) 

As the spikes damp out the solution finally enters the small-signal phase 

Vu - ¥ m -» 0. (55) 

The small-signal solution is well known, ' 6,13 ' 2 9 ' 31 so there is no need 
to discuss it here, but we give it for ready reference: 

r?(r) = 1 + Ae~ hu * T [n cos (ilr + tp) - ±co£ sin (Ut + <p)] 

p(r) = /3[1 + Ac~ }MT sin (At + *)] (56) 

(f = tap - i(^) 2 

where A is an arbitrary real small amplitude and <p is an arbitrary real 
phase. Exactly the same small-signal solution is obtained from (25), 
thereby justifying the approximation (22). It is now easy to see that 
when (16) breaks down the frequency ft of the small-signal solution be- 
comes imaginary and there are no oscillations. This may be the case in 
the gas laser, where we may have £ ~ 1, «. 1, and co > 4/3. 

IV. A NUMERICAL EXAMPLE 

Before attempting to apply our formulas to the analysis of spiking 
patterns we wish to discuss their accuracy with the aid of a machine 
calculation. For numerical integration the rate equations (11), (12) are 
best written in the form 

± = «(£ — x — xe") 

(57) 

y = ac " + x - 1, 
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where 

X = -q 

(58) 
y = In p. 

As already pointed out following (17), the initial condition on x,y is not 
critical; it is convenient to start the solution on the adiabatic solution 
(17), so we take 

?/(0) = hi a, x(0) = 0. (59) 

The constants will have the values 

to = 7.12 X 10~ 6 

€ = 5 (60) 

<7 = 2 X 10~ 9 

which are typical for a ruby laser. The numerical integration of (57) 
has been performed on the IBM 7090 computer using Hamming's 32 
predictor-corrector method. The program provides for automatic halving 
and doubling of the integration interval under the control of the frac- 
tional error of the increment and a preselected tolerance (5 X 10" ). 
The stability and accuracy of the program for this problem were checked 
using initial conditions which lead to the small-signal solution (56). 

The results are shown in Fig. 3. The solid curve and scale on the right 
give y(r), while the dotted curve and scale on the left give x(t) for t 
in the range 3 X 10 4 to 4 X 10 4 . The origin for r is completely arbitrary. 
From (26), (58) and (60) we have 

2/ = * + ln/3 = *-f- 1.387. (61) 

We are concerned primarily with y, since the inversion is not ordinarily 
observed experimentally. The main features of the computed results are 
summarized in Table I, which lists the values of t„ , y n , and m n or M n 
for n = 1, • • ■ , 9 for the first nine extrema in the notation of (34). These 
values were obtained from the computed points by fitting a parabola to 
the three points nearest the extremum. The spacing of the computed 
points was sufficiently small (At = 0.005 X 10 4 ) that in all cases all 
three points lay well within the validity of the parabolic approximations 
(38) or (41). 

We shall now attempt to calculate the information of Table I by the 
formulas of the preceding section. In every case we shall indicate the 
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Fig. 3 — A machine calculation of y(r), solid curve with scale on the right, and 
x(t), dotted curve with scale on the left, satisfying the rate equations (57) for 
<i> = 7.12 X 10 -6 , £ = 5, a = 2 X 10 -9 . The important results are summarized in 
Table I. Here x(t) represents the inversion and y(r) the logarithmic light output. 



Table I — Computed Results of Numerical Example 

Summary of results of machine solution to the rate equations for w = 7.12 X 
10-°, £ - 5, a = 2 X 10" 9 . Also obtained were x max - 1 = 0.02830, y max = 
0.0282. 



n 


Tn 


Vn 


tn n 


M n 


1 


3.1344 X 10* 


-14.570 


— 




2 


3.2521 


4.240 




0.01416 X 10 4 


3 


3.3642 


-12.815 


0.0532 




4 


3.4806 


4.179 




0.01415 


5 


3.5906 


-11.795 


0.0529 
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3.7013 


4.053 




0.0151 
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3.8100 


-10.986 


0.0531 
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3.9186 


4.064 




0.0152 
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4.0239 


-10.303 


0.0531 
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correct value as computed by machine on the right in parentheses. Ac- 
cording to (19) 

(02) 
y x = Inpi = -14.09 (-14.57). 

Thus we overestimate iji by 0.48 or 3.3 per cent. From (61 ) the theoretical 
value for \^i is 

*i = -15.48. (63) 

Putting ¥« = ^i and ^r u = ^2 in (30) gives 

* 2 = 2.91, (64) 

or 

;/ 2 = 4.30 (4.240). (65) 

Thus we overestimate y-> by 0.06 or 1.4 per cent. For the excursion of y 
we obtain 

,/, - , h = *., - *, = 18.39 (18.81), (66) 

which is an underestimate by 0.42 or 2.2 per cent. Henceforth we need 
not consider ^, but only y. From (50) using ^.„ — ^„. = y-> — 2/1 we 
obtain 

2/3-2/1 = 0.965 (1.355) (67) 

7/3 = -13.12 (-12.815). (68) 

We are underestimating y A by 0.30 or 2.2 per cent. From (52) using 
yfr M — <& m = 7/2 — 2/3 we obtain 

y 2 - y 4 = 0.0541 (0.061) (69) 

2/4 = 4.25 (4.179). (70) 

By repeating steps (67) and (69) we could obtain values for any number 
of succeeding maxima and minima 

?/ 5 - y-s = 0.881 
7/5 = -12.24 

y 4 - y/6 = 0.0523 

2/6 = 4.19 

We observe from Table I that 2/4 - 2/o is unusually large (0.126), 



(1.020) 


(71) 


(-11.80) 


(72) 


(0.126) 


(73) 


(4.053). 


(74) 
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while y 6 — y s is negative ( —0.01 1 ) ; the average of these is 5(2/4 — y&) = 
(0.058) in better agreement with the machine-computed y 2 — y\ and 
with theory. Thus we suspect that the discrepancy in (73) is not sig- 
nificant. There is apparently a slight inaccuracy in the computed solu- 
tion on the third and fourth spike. One might expect that noise of 
various kinds would have a similar effect on real lasers, making the damp- 
ing from spike to spike unreliable. Averaging over several spikes, how- 
ever, would still give a meaningful value for Ay M , as it does in our com- 
puted spike pattern. 

The durations of all the minima should be given by (39) 

m = 0.0520 X 10 4 (0.0531 X 10 4 ). (75) 

We see in Table I that m is constant within 0.4 per cent and the average 
over the four minima is 0.0531 X 10 4 . The durations of the maxima will 
be calculated from (42), using for Vm — &m the values of the excursions 
just calculated, with the minimum preceding the maximum 

Mi = 0.0127 X 10 4 (0.01416 X 10 4 ) 

Mt = 0.0130 (0.01415) (76) 

M 6 = 0.0134 (0.0152). 

Here we have underestimated the spike durations by 11-12 per cent. 
The intervals are given by (46) 

h. 4 = 0.226 X 10 4 (0.2285 X 10 4 ) 

7 4 . 6 = 0.220 X 10* (0.2207 X 10 4 ). 

The agreement here may be considered perfect. According to (53) we 
should have 

by*/ 1 - 2.37 X 10 -5 (78) 

for all consecutive spikes; the computed results are as follows: 



M(4 



2-*4 




(2.66 X 10 -5 ) 


4->6 




(5.70 X 10 -5 ) 


0->8 




(0.51 X 10~ 5 ) 


6) + (6- 


-8)] 


(2.63 X 10" 5 ). 



(79) 



The last line is considered more significant than the discrepancies in 
the third and fourth lines. Thus our theory underestimates Ay^/I by 
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11 per cent. Finally, from (35) 

fc»a, = -W - 1 = 0.0288 
(0.0282) (0.0283), 

which shows that the approximation (21) is excellent. 

We may conclude from this confrontation between our formulas and 
a machine computation that the accuracy is probably adequate for the 
analysis of experimental spike patterns. 

V. APPLICATION OF THE FORMULAS 

Experimental data on spiking are ordinarily obtained in the form of an 
oscilloscope trace proportional to the light output. The trace is propor- 
tional to p(t), the number of photons in the laser mode of the cavity, 
but it is usually considered impractical to calibrate the equipment so 
as to obtain the absolute value of p(t). And even if p(t) could be meas- 
ured absolutely it would not fix the absolute value of p(r) defined in 
(9). We shall assume that only relative values of p(r) can be measured. 
The most significant measurements of p are the peak-to-valley ratios 

R = PM/ P ,n, (81) 

the ratio of the maximum light output of a spike to the minimum output 
in a neighboring minimum. In analyzing a spike pattern with many 
spikes, some convention must be adopted on which minimum to choose. 
From (26) and (61) 

*m - *« = Vu - Um = In R. (82) 

Thus the excursions of the model particle in the potential V(¥) are 
uniquely fixed by the data. It follows from (30) that both ¥« and V m 
are fixed by R 

*,, = ln(lnfl) (g3) 

¥ m = In (In R) - In R. 

Even M>i , the hypothetical minimum not actually observed, can be de- 
termined by extrapolating the other minima ty m . Thus all of the extrema 
may be regarded as immediately fixed by the data. From the extrema 
alone we obtain 2 relations satisfied by the three parameters w,£,o-, 
namely (19) and (50) or (52). 

It is now clear that a single spike pattern does not contain enough 
information to determine w,£,<r. The reason for this is that we cannot 
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make use of the measured durations or intervals, since we do not know 
the ratio t p between real time / and dimensionless time r. We might 
think that we could obtain information from time ratios such as (M/m) 
or (m/I). It turns out that we can, but it is a different kind of informa- 
tion, the kind we have called "qualitative" in Section I. Consider the 
three quantities 

.-1 = (M/I) Infl (84) 

B = (m/7)[(ln /? - l) 5 + (In R - 1 )" J ] (85) 

C - (Atftf/A¥ m )(ln# - 1), (86) 

all of which can be determined immediately from the data since they 
depend only on ratios. It is to be expected that C might have to be 
averaged over several spikes to get a meaningful value. Except for this 
difficulty, values of .1, B, and C may be calculated for every spike in 
the data after adopting some convention to handle the usual ambiguity 
in the definitions (84), (85), (80). One such convention is illustrated 
by the following example : 

In R-2 = * 2 - * 3 

A, = (MJI M ) In/?.. 

(87) 
B 2 = (w ;i //,_ 4 )[(ln R 2 - 1 Y + (In /?., - 1 P'] 

('■2 = [(*j - *4>/(*i - *i)](In ff 2 ~ 1). 

From (42), (40), (50) and (52) we find the simple relations 

A = B = C = I (88) 

which do not involve the parameters a>,£,o\ Therefore (88) should hold 
even if w,£ are slowly-varying functions of time within our adiabatic 
approximation. We call the relations (88) consistency relations, since 
they can be used to determine whether data are consistent with the 
rate equations. If (88) is satisfied reasonably well, at least in an average 
sense over the spikes, it is a foregone conclusion that a reasonable fit to 
the data can be obtained from the rate equations. The converse is also 
true: if (88) is not satisfied the data cannot be fitted from the rate 
equations. 

The importance of observing the minima (valleys) as well as the 
maxima (peaks) in a spike pattern is abundantly clear. Without the 
valleys we cannot determine R, m, or A^,„ , all of which appear in the 
consistency relations. Therefore a great deal of information is lost unless 
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the valleys can be seen above noise. This presents some difficulty, 
because ordinarily R » 1 is much too large to be measured from an 
oscilloscope trace that responds linearly to the light output and contains 
both the peaks and the valleys. In spike patterns the valleys usually 
just correspond to the noise level in the experiment. It is not our purpose 
to go into experimental details except to point out that with care the 
valleys should be observable. The basic experimental requirement is 
that the acceptance cone of the detector should correspond to the radia- 
tion cone of the laser so as to exclude extraneous light from the pump 
and the spontaneous emission of the laser medium. If it is not possible 
to measure the valleys, there is still one consistency relation that can 
be applied to the peaks alone. From (46), (82) and (83) we have 

m R = i/M 

-*„, = (I/M) - In (I/M). 

The use of (46) is equivalent to assuming A = 1. Thus the minima can 
be calculated from I/M if we assume A = 1. We could obtain the 
durations m of the minima by putting 5=1. The damping of the 
minima would not be given very reliably by (89), but can be obtained 
from (86) by putting C = 1. Thus we can deduce R, m, and A^ m from 
the peaks alone, but we lose all of our consistency relations (88). How- 
ever, from (89), (42), (39) and (35) it follows that 

D= 1, (90) 

where 

D = (MviW4)[(//M) " if/KUM) - In (I/M) - \f (91) 

can be determined from the peaks if the time resolution is good enough 
to give a good value for (My max ). If (90) is satisfied, it is probably good 
evidence that the laser obeys the rate equations, and a rate equation 
analysis is meaningful. However, if all that is desired is to apply the 
rate equations blindly to obtain quantitative information, it is not neces- 
sary to measure M^ max . All of the quantitative information in a con- 
sistent spike pattern can be deduced from I, M, and A^*, . 

We now consider practical ways of obtaining quantitative information 
from spike patterns. The one-pattern method is to measure t p by an inde- 
pendent experiment. The measurement of the cavity losses has been 
discussed by several authors. 26,33 Suffice it to say here that t p can be 
measured from the dependence of the threshold flash energy for produc- 
ing laser action on the temperature and on losses deliberately introduced 
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into the laser mode. Once t p is measured all of the times m, M and I 
become known in dimensionless time. From (39) 

up = w (£ - 1) = 8/m 2 , (92) 

and from (53) 

co£ = \&* M /L (93) 

These equations can be solved for co,£ 

£ = ym/(ym — 1) 

co = 8(ym — lj/ffl, 

where 

7 = A(m/J)A*„ (95) 

is independent of t p . Now <j can be obtained from ( 19) and (89 ) 

a = ±p(a>P)\l/M)e- ( " u) (96) 

with M = Mi and / = 7 2 -4 • This procedure makes use of the duration 
m of the valleys but not the peak-to-valley ratio R. If only the peaks 
are observed we write instead of (92) 

up = „(g - 1) - (8/MI)/[l - (M/I)]. (97) 

Solving (93) and (97) for co,£ gives 

£ = SM/(8M - 1) 

(98) 
co = 8(5.1/ - l)/(/ - M)M, 

where 

8 = A[l - (itf//)]A¥* (99) 

is independent of /„ . This kind of analysis can be applied to every spike 
in a spike pattern. It may be found that £ and especially co vary slowly 
from spike to spike, which is permissible within our adiabatic approxima- 
tion. If the data satisfy the consistency relations (88), the same values 
of co,£ will be obtained from (94) and (98). The greatest weakness of 
this method is that (94) fails completely if ym ^ 1 ; (98) fails if 8M ^ 1 . 
Thus a great deal depends on the accuracy of the formulas as well as 
that of the measurements of t p and y or 8. It follows that the method will 
fail whenever £ » 1. 

As an example we shall apply this method to the machine calculation 
considered in Section IV. We consider the solid curve of Fig. 3 to be the 
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data, which implies that t p is known and both the peaks and valleys 
have been studied. The relevant numbers have already been given in 
parentheses in (69), (75), (76) and (77); we repeat them here without 
parentheses 

A* a, = i j2 - y 4 = 0.061 

m = 0.0526 X 10 4 

M = M 2 = 0.0142 X 10 4 

I = / 2 . 4 = 0.2285 X 10 4 . 

From (94), (95) and (96) we obtain 

k = 3.4 (5.0) 

co = 12 X 10" 6 (7.12 X 10~ 6 ) (101) 

a = 10 X 10~ 9 (2 X 10~ 9 ). 

The lack of accuracy is primarily due to the fact that £ = 5 is a little 
too large to give good results with this method. Using only the peaks, 
we obtain from (96), (98) and (99) 

£ = 2.9 (5.0) 

co = 14 X 10 -6 (7.12 X 10~ 6 ) (102) 

a = 8 X 10" 9 (2 X 10~ 9 ). 

We now describe a two-pattern method which does not require the 
measurement of t v . In fact, it yields a value for t p and may in some cases 
give better results for co,£ than the one just described. It is based upon 
observing the valleys in two or more spike patterns for which the 
relative values of co and £ are known. We must assume that the linewidth 
A v is known as a function of temperature. Let us suppose that we meas- 
ure the valley durations m and m' in two spike patterns in which the 
ratios (£'/£) and (co'/co) are known. We know (£'/£) from the tempera- 
tures at which the patterns were obtained. We can obtain co'/co from the 
relative pump intensities at the times when spiking occurred. It im- 
mediately follows from (92) that 

_ (cow /com ) - 1 , s 

* (com 2 /"'™' 2 ) - (€70 • ; 

This result is meaningful providing that the data give (cora 2 /co'ra' 2 ) lying 
between (£'/£) and unity. Once £ is determined co can be calculated from 
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(92) and (93) 



9 U- 

a) = 



1) (!lj (A *„) 2 , (104) 



32 e 

and a is again obtained from (96). 

As an example of this method we shall apply it to two machine-calcu- 
lated spike patterns, one of which is that of Fig. 3 with the parameters 
(60), and the other has the parameters 



I' - 4 (105) 

a = a. 

From the machine-calculated pattern we find 

m ' = m3 ' = 0.0614 X 10 4 . (106) 

The value of m = m 3 is given in (75). We assume that the ratio 

$'/£ = 0.8 (107) 

is known from the temperatures at which the data were taken, and the 
ratio 

(co?n7a/m /2 ) = 0.748 (108) 

is calculable from (oj/V) and the two valleys. We now obtain from (103) 
the value 

$ = 4.85 (5.0). (109) 

From (104) we now obtain 

w = 9.0 X 10~ 6 (7.12 X 10" 6 ) (110) 

using the values in parentheses from (75), (77) and (69). From (39), 
(109) and (110) we obtain the absolute value of m 

m = 0.048 X 10 4 (0.0531 X 10 4 ). (Ill) 

It follows that t p is given by 

t P = m exp /m, (112) 

where m exp is the measured duration in laboratory time. We conclude 
that the two-pattern method is to be preferred to the one-pattern method 
as a general approach to spike analysis. It should be mentioned that 
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there seems to be no two-pattern method involving only the peaks which 
is sufficiently accurate to give meaningful results. 

The most convenient experiments to perform involve changing the 
pump intensity l/t„ while all other parameters remain fixed. If t„ « U , 
as is usually the case in flash-pumped lasers, £ will be independent of t g 
and u) will vary as \/t„ . Varying w does not give quantitative informa- 
tion such as we obtained from varying £ in the two-pattern method. 
Nevertheless, it is of interest to consider what effects are to be expected. 
From (30), (26) and (19) we have roughly 

p.,~-w„, (114) 

~/3 In C/3-\/co/3/2cr). 

Thus the maxima in p depend only logarithmically on pumping power. 
It follows from (9) that the observed light output proportional to p 
should vary as 

V « 1/V (115) 

The time intervals m, M, and I should vary as 

m oc M °c / « t}. (116) 

VI. SUMMARY 

We have now outlined a comprehensive program for the study of lasers 
by means of their spiking patterns. The rate equations have been formu- 
lated in terms of the light output and the atomic inversion and five 
physical parameters in (1) and (2). We have written the equations in a 
general form valid for three- and four-level systems and taking into ac- 
count statistical weights and thermalization of the upper laser level. 
These equations were then reduced to dimensionless form in terms of 
three parameters in (11) and (12). All of the properties of the spiking 
solutions of experimental interest were then deduced analytically by 
means of the mechanical analogy shown in Fig. 2, in which a particle 
moves in a potential well in a viscous medium. 

The formulas obtained were illustrated and tested for accuracy against 
a machine-computed solution to the rate equations. The value of the 
formulas was confirmed by this comparison, and we went on to discuss 
their application to experimentally observed spiking patterns. Four rela- 
tions were given whereby spike patterns can be tested for consistency 
with the rate equations. These consistency relations do not contain any 
parameters, only ratios of times and of light outputs. 
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Two methods were described for obtaining quantitative information. 
In the one-pattern method everything is deduced from a single spike 
pattern, but it is necessary to measure the photon lifetime l p in the cavity 
by independent experiments. It is possible to apply this method to data 
in which only the peaks are observed. It is emphasized, however, that 
the observation of the valleys in light output should be perfectly feasible 
and very much worthwhile. In the two-pattern method all the parameters 
and l p are deduced from two patterns obtained at different temperatures. 
In this method, the more accurate of the two, it is essential that the val- 
leys be observed. 

Our objective has been to provide the researcher with a set of tools for 
applying the rate equations to experimental spiking data. The spiking 
phenomenon in solid state lasers can be utilized in research now that it is 
beginning to come under good experimental control. We have tried here 
to point out what kind of spiking data is needed, and what additional 
information is needed, to get the maximum information from spiking. 
Our analysis applies to the sharp spike region of time in which the rate 
equations are highly nonlinear. It complements the well known small- 
signal analysis in which the rate equations become linear. The entire dis- 
cussion is based upon the rate equations of Statz and De Mars, which we 
regard as reasonable, relevant and widely accepted. We have not gone 
into the derivation of these equations or the modifications that have been 
proposed or might be proposed. We prefer to leave that to the future, 
when we may reasonably expect that systematic spiking studies will 
have clearly revealed the adequacy or inadequacy of these equations, 
and indicated the direction which these modifications must take. 
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